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^ : Abstract 
l> . 

We present exact, explicit, convergent periodic-orbit expansions for individual 

■ energy levels of regular quantum graphs. One simple application is the energy 

■ levels of a particle in a piecewise constant potential. Since the classical ray 

trajectories (including ray splitting) in such systems are strongly chaotic, this 

result provides the first explicit quantization of a classically chaotic system. 
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I. INTRODUCTION 



Within the framework of semiclassical periodic-orbit theory the quantization procedures 
for integrable and chaotic systems differ substantially. An integrable system may be quan- 
tized using EBK theory [0]. The set of integrals 

Pf,dq^ = h{ni + fii), i = l,...,N, (1) 

J Ci 

extended along the fundamental cycles Ci of the A^-dimensional phase-space tori yield 
the (semiclassical) quantization conditions for every action variable Jj. Here the rij's are 
integer quantum numbers and the /ij's are Maslov indices. Although not exact in general, 
the quantization condition does (implicitly) produce individual energy levels -Eni,...,njv 
that can be labelled, one by one, with the A^ quantum numbers ni, . . . , n^r. This procedure 
differs markedly from the chaotic case where the focus is not on individual energy levels but 
on global characteristics of the spectrum. For instance, instead of finding individual energy 
levels as in (|ID, periodic-orbit quantization schemes for chaotic systems, such as Gutzwiller's 
trace formula compute the density of states 

oo 

p{k)=Y.Hk~k,), (2) 

from which individual energy levels are extracted indirectly as the singularities of p. In a 
chaotic system the only available classical input are the periodic orbits of the system and 
the density of states (0) is computed according to 

1 oo 

p{k) ^ -p{k) + - Im^T,(E) Y: A;{E) e-^^(^). (3) 

Here p(fc) is the average density of states, Sp{E), Tp{E) and Ap{E) are correspondingly 
the action, the period and the weight factor of the prime periodic orbit labeled by p, and 
u is the repetition index. Again, the scheme (^ is not usually exact. More seriously, 
however, in contrast to (|1]) it fails to produce individual energy levels in the form = ...". 
The difference between (|I]) and (^ cannot be emphasized enough. While (|ip allows us to 
"pick and choose" a particular energy eigenvalue, in the chaotic case all of the eigenvalues 

2 



have to be computed according to (|^), and only a subsequent non-analytic inspection and 
counting procedure allows us to focus on an individual energy level. There is, however, a 
class of quantum chaotic systems, regular quantum graphs 0, which are explicitly solvable 
analytically 0, i.e. exact periodic-orbit expansions of the form = ..." exist. The 
purpose of this paper is to expand considerably with respect to the work presented in 
and to present a thorough discussion of our methods and their validity. 

The organization of this paper is as follows. In Sect. II we extend the theory of quantum 
graphs [0-0 to include dressed graphs, i.e. quantum graphs with arbitrary potentials on 
their bonds. In Sect. Ill we define regular quantum graphs and present explicit, convergent 
periodic-orbit expansions of individual eigenvalues. These expansions are not just formal 
identities; the periodic-orbit expansions presented in Sect. Ill converge, and converge to the 
correct eigenvalues. In Sect. IV we present a worked example of a simple quantum graph 
whose spectrum is computed in three different ways: (i) numerically exactly, (ii) via the 
explicit periodic-orbit expansions presented in Sect. Ill and (iii) via numerical integration 
using an exact trace formula for the density of states. The results of the three methods 
agree. This proves validity and convergence of our approach. In Sect. V we summarize our 
results and conclude the paper. 



II. DRESSED QUANTUM GRAPHS 

A quantum graph consists of a quantum particle whose motion is confined to a one- 
dimensional network of Nb bonds Bij connecting Ny vertices Vi. An example of a graph 
with six vertices and ten bonds is shown in Fig. 1. The topology of a given graph is fully 
characterized by its connectivity matrix Cij: 

{1, if Vi and Vj are connected, 
(4) 
0, if they are not. 

Every bond Bij which connects the vertices Vi and Vj, supports a solution ipij{x) of the 
Schrodinger equation 
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\lx " ^'^J = Ei)ij{x). (5) 

Here < x < Lij is the coordinate along Bij measured from Vi to Vj, and Lij = Lji is 
the length of the bond. A constant, real, skew symmetric matrix Aij = —Aji, which plays 
the role of a magnetic field vector potential, is sometimes introduced as a tool for braking 
the time-reversal symmetry, which, in turn, is known to affect the statistics of the level 
distribution ||,||. 

In this paper we generalize the Schrodinger operator in (||) by adding potentials Uij{x, E) 
on the graph bonds. We call this generalization "dressing of the graph" . While in general the 
potentials Uij{x, E) may depend on the bond coordinate x and the energy E in an arbitrary 
way, we restrict ourselves here to the scaling case 

Uij{E) = XijE, Xij = Xji, (6) 

which allows us to introduce physical parallels between quantum graphs and ray-splitting 
systems [^U]-[T^. A quantum graph with the potentials (|^) on its bonds can also be viewed 
as a generalized step potential, such as the one shown in Fig. 2 (a). These potentials were 
studied before in great detail in connection with Anderson localization [|l^. Potentials of 
this type can be represented by a linear graph, such as the one shown in Fig. 2 (b). Scaled 
potentials such as @ cast the Schrodinger equation into the form 

dx 



-'^-^^A M^) = 0tJE^p^,{x), (7) 



where the parameters Pfj = 1 — Ay, Pij = (3ji are defined on the corresponding bonds Bij. 

Depending on whether the energy E = k"^ oi the particle is above or below the scaled 
potential height Uij{E), the solution of equation (|^) on the bond B^j is either a combination 
of the free waves, 

f^ij k-\-A.ij')x ^i(^f3ij k-\-A.ij )x 

^ij{x) = aij j= — + bij , Xij < 1, (8) 



or a combination of the tunneling solutions. 



ipij{x) = aijc^-^^i^+'^^i^'' + %e('''^'=+'^>^>, Xij > 1, (9) 

where the factors {Pijk)~^^'^ in the propagating waves (H) were introduced to insure proper 
flux normahzation [|14| . Due to the scahng assumption, there is no transition between these 
two cases as a function of E. From now on we shall assume that the energy E is kept above 
the maximal scaled potential height, 

Ay<l, i,j = l,...,Nv, a, ^0, (10) 

which will allow us to exclude the tunneling solutions (^). At every vertex Vi, the bond wave 
functions satisfy the continuity conditions 

\x=o= 'fiCij, i,j = 1, ...,Nv (11) 

and the current conservation conditions 
( d \ 

5Z y-i-^ - Ajj i'ijix) \x=o= -iKfi, hi = 1, Ny. (12) 

Here ipi is the value of the wave function at the vertex Vi, and the Aj's are free parameters 
of the problem for which the scaling is introduced as 

Xi = A°A;. (13) 



The conditions (|TT]) and (|T2D are consistent only for a discrete set of wave numbers fc„, which 
defines the spectrum of the quantum graph problem. Since ipij{x) and ipji{y) represent the 
same wave function on the bond connecting the vertices Vi and Vj (the only difference is 
that X is measured from vertex Vi and y is measured from vertex Vj) we have 

ijjji{Lij - x) = tpij{x). (14) 

Using (H) we obtain 

^i{-/3ijk+Aji){Lij-x) ^i{l3ijk+Aji){Lij-x) 
Ipji {Lij -x) = aji j= h bji j= = Ipijix). (15) 

Therefore the coefficients aij and bij are related according to 



a^.. = 6.^.e'(^»i^+^»^)^'^ bji = aije'^-'^^^'^+^^i^^'J . (16) 

The coefficients aij and aji, {pij and hji) are considered different. This imphes that the bonds 
of the graph are directed. Equations ( [T6| ) can be written in matrix form, 

a = P]D{k)h, (17) 

where a and h are the 2A'^B-dimensional vectors of coefficients, D is a diagonal matrix in the 
2Nb X 2Nb space of directed bonds and 

P=[ h (18) 

Vn, y 

where 1 Nb is the A'^g-dimensional unit matrix. Exphcitly we have 

A,,p,(A;) = M.<?e*^''-'+^-)''-. (19) 

The pairs of indices (ij), (pg), identifying the bonds of the graph F, play the role of the 
indices of the matrix D{k). Alternatively the wave function can be written as a linear 
combination of plane waves scattering off Vi. An incoming wave with normalized flux on the 
bond Bj/i gives rise to a partial wave contribution scattering into bond Bij according to 



' Hij 1^ y Pij 1^ 



Here crjj'(fc) is the matrix element of the vertex scattering matrix a^'^\k), which distributes 
the incoming flux on bond Bj/i into the bond Bij. The wave function ipijlxj) on the bond 
Bij is a superposition of the partial waves (|20D with amplitudes a^j/ corresponding to the 
incoming flux on the bond Bj/i towards the vertex Vi, i.e. 

r 

Using the representation (H) of ifjij in (^) and comparing coefficients yields 

% = E^S'«^.'- (22) 



Substituting (^) into the boundary conditions, we obtain the vertex scattering matrix 



-fr - 4!..' = I + 1 C^^^^^'' (23) 



with 

t;. = EA,a,. (24) 

i=i 

We see that in the scahng case the matrix elements aj*]/ of the vertex scattering matrix 
cr'^*) are ^-independent constants. The matrix element crj*] has the meaning of the reflection 
coefficient from the vertex Vi along the bond -Bj^, and the elements crj*j,, j ^ j' are the trans- 
mission coefficients for transitions between different bonds. Equation p^ ) can be written 



as 

b = fa, (25) 

where 

= Tij^nm ^inCjiCnm'~^j[m' (26) 

Equations (|17|) and (pSf) together result in 

a = S{k)a, (27) 

where S{k) (the total graph scattering matrix) is given by 

S{k) = D{k)T (28) 

and D = PDP and T = PT. The consistency of the system of linear equations ( p?!) requires 
the spectral equation 

A(A;) = det [1 - S{k)] = (29) 

to be satisfied. This condition defines the set of allowed momenta 

The density of the momentum states of the dressed quantum graph is given by 



p{k) = Y.5{k-kn), (30) 

n=l 

where the A;„'s are the solutions of (^91) . An exact periodic-orbit expansion for p{k) can be 
obtained directly from the spectral equation (pOf) as follows The logarithmic deriva- 

tive of (|29D is singular at every one of its roots. Between roots, the phase of the spectral 
determinant varies slowly such that 

p(k) = p(k) - - lim Im 4r In det [1 - ^ (A; + ie)] . (31) 

Using the well-known identity 

oo -I 

In det[l -S] = -TrJ2 - S", (32) 



n=l ^ 

we obtain 

pik) = pik) + - hm Im 4t y - Tr \S{k + «e)]" . (33) 

Since the matrix indices of ( pSD correspond to the vertices of the graph, the trace of the 
nth power of the scattering matrix can be interpreted as a sum over all closed connected 
sequences consisting of n bonds 0-0] • Classically, these periodic connected sequences of n 
bonds Bij correspond to the periodic orbits traced by a point particle moving on the graph. 
Geometry and proliferation properties of the periodic orbits are determined completely by 
the topology of the graph. 

The behavior of the periodic orbits on graphs exhibits the typical features of chaotic 
systems. The meaning of classical chaoticity on graphs is well defined, as demonstrated in 
the following. A classical graph system consists of a graph F and a point particle moving 
along its bonds, which scatters elastically at every vertex Vi along the direction of any of the 
bonds emanating from this particular vertex, with different probabilities. The probability 
amplitudes for every scattering channel can be obtained in the short wavelength limit from 
the quantum mechanical transition amplitudes defined at every vertex Vi by the correspond- 
ing scattering matrix crj'j,. In the scaling case, the matrix elements crj*j/ are fc-independent 



constants and thus do not depend on h at all. Therefore the same matrix elements determine 
both the quantum and the classical scattering probabilities. 

For every given graph F the global average rate of exponential proliferation of periodic 
orbits, the topological entropy A, is given by 

A=lin,!!iMM, (34) 

Z^oo I 

where / characterizes the length of the periodic orbits (for instance their code lengths) and 
is the total number of periodic orbits of length < /. The number of possible periodic 
orbits increases exponentially with their lengths, (or, equivalently, the number of scattering 
events) with a rate A which depends only on the topology of the graph. Since the phase 
space of the system is bounded, the dynamics of such a particle is mixing [Q. 

Since we are focusing on the case Xij < 1, most of the classical periodic orbits on a 
graph are above-barrier reflection orbits as illustrated in Figs. 2 and 3. In the context of ray 



splitting these orbits are also known as non-Newtonian orbits [15-17|. The inclusion of all 
non-Newtonian orbits in our periodic-orbit expansions of individual eigenvalues discussed 
below is crucial for rendering these expansions exact. 
Traversing the bond Bij contributes the amount 



Sij = / iPijk + Aij) dx (35) 



to the total action of the trajectory traced by the particle. These actions appear in the phases 
of the exact wave functions (H). This connection means that the semiclassical (eikonal) form 
is exact for the quantum graph wave functions. More importantly, the amplitudes e*'^*-'' 
determine the matrix D{k), and hence the scattering matrix S{k). As a consequence, the 
"closed bond sequence expansion" (|33|) can be written explicitly as a periodic-orbit expansion 
in terms of the phases (p5|), 

1 °° _ 

p{k) = p{k) + - Rej:T,{k) E A;e^^'^^''\ (36) 



where Sp is the action of the prime periodic orbit p composed of the partial actions Sij of 



equation (|35|) accumulated along the periodic orbit p, and Tp{k) = dSp{k)/dk. The first term 



in (^6D describes the average behavior of the density of states while the second represents 
the fluctuations around the average. The amplitude of every periodic orbit p contains the 
constant factor exp{iJ2ij^ijLij)- This factor can be absorbed into the weight factor Ap. 
Thus, defining the reduced classical bond actions 

Sij = PijLij (37) 

and the total reduced action Sp accumulated along the periodic orbit p, 

S',= E si (38) 

ij along p 

the final periodic-orbit expansion for the density of states for scaling systems can be written 

as 

-1 oo 

p{k) = p{k) + - Re E E a; (39) 

^ p u=l 

In contrast with the expression (|39| ) for the density of states is exact; the action lengths 
Sp and the weight factors Ap are /c-independent constants. 
The staircase function 

oo 

N{k) = E - kn) (40) 

n=l 

is obtained by direct integration of (0). Using (p9D , N{k) can be expanded as 

N{k) = N{k) lim Im In det [1 - S{k + ie)] = N{k) + - Im E E — e*^^°^ (41) 



Just like (|39|) this expansion is exact. The first term represents the average behavior of the 
staircase; the second term describes the fluctuations around the average. 

III. REGULAR QUANTUM GRAPHS AND EXPLICIT SPECTRAL FORMULA 



Since the scattering matrix ( p8D is a unitary matrix, its eigenvalues have the form S[ 
Qi9i{k)_ Therefore the spectral determinant (EUf) can be written as 
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2Nb 

A{k) = n [i - e'^'^''^ 
1=1 

— 2e*®o(fc) 



2Nb 



1_ ^e^e.w^^^giEil? 



1=1 



Nc-l 

cos[eo(A;)]+ £ (-iycos[e,(A;)] 



(42) 



where 



2Nb 



eo 



1=1 



is the total phase of the spectral determinant, the 6j's in 



(43) 



are linear combinations 



(sums and differences) of the phases 9i{k) and Nc = 4^s/2. Evaluated directly, the spectral 



determinant is a polynomial of the (complex) matrix elements (19) with coefficients that 
are determined by the matrix elements (|26|) . Factoring out the total phase (^31) of this 
polynomial, we obtain the spectral equation in the form 



cos 



(44) 



where 



Nr 



^{k) = E cos(r2jA; — Tfji 



(45) 



i=l 



Here, based on the reduced bond actions defined in (^), Sq = Y,ij CijS'^^ is the total reduced 
action length of the graph, the frequencies Qi < Sq are sums and differences of the reduced 
bond actions 5"°- and 70, 7i are constants. For a general graph F it is difficult to calculate 
the precise number of cos-terms A^r in (|45|) . But an upper limit is given by the number of 
possible linear combinations of the A''^ reduced bond actions Sfj. Since there are 2^ 
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ways of picking j actions out of Nb possible ones and combining them with "+" and "— " 
signs, we obtain 

'Nb\ Nb / Nb' 



Nb 

iVr < E 2'' 



Nb 

-1 + E 



(46) 



A graph F is called regular, if the condition 

a <l 



Nr 

E 

i=l 



(47) 
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is fulfilled. In case the condition (|^) is satisfied, the spectral equation ( ^41) can be immedi- 
ately solved to yield the following implicit equation for the eigenvalues 

^ I ( arccos[$(fc„)], for n + /i even 

kn = ^[n + fJ. + ^o] + ^\ (48) 
"-^0 *-^o TT — arccos[<P(/c„jJ, tor n + yU odd, 



where /i is a fixed integer, chosen such that ki is the first positive solution of (pf). Equation 
(^) implies the existence of separating points 

kn = ^{n + fi + -fo + l) (49) 



in the spectrum fc„ of (0). Because of (]47|) the points A;„ are never solutions of (^81) . They 
act as separators between fc„ and Since the second term in ( ^Sf ) is bounded by h/Sq, 

the deviation | A;„ — A;„ | never exceeds vt/S'q for any n. We emphasize, that the separators 
kn do not coincide with the average values A;„ of the roots A;„. The explicit decomposition of 
the roots fc„ into an average part A;„ and a fluctuating part A;„, A;„ = A;„ + fc„, can be obtained 
from the following equivalent formulation of (^HD, 



~ 5*0 



1 

^ + + 7o + 2 



CO 
On 



|arccos[$(A;J] - II . (50) 



'0 

This form of kn together with the boundedness of the second term in ( ^0|) proves rigorously 
that N{k), p{k) are of the form 

iV(A:) = ^A: + iV(0), p{k) = ^ = ^- (51) 

TT ak TT 

Since $(A;) contains only frequencies smaller than Sq, every open interval J„, = (A;„_i,A;„) 
contains one and only one root, i.e. kn- Moreover, if (^) is fulfilled, the allowed zones 
Zn C In where the roots kn can be found narrows to 

(71 71 \ 

-^(n + /i + 7o + u),-^(n + /i + 7o + l-u) 1 , (52) 

where u = arccos(a)/5'o. Correspondingly, there exist forbidden regions Rn, 

(71 71 \ 

(n + /X + 7o + 1 - m) , (n + /X + 7o + 1 + m) 1 , (53) 
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where roots of ( ^8]) can never be found. Note that A;„ G -R„. In the hmit a — > 1 (-u — » 0), 
the width of the forbidden region i?„ shrinks to zero, and the allowed zone Zn occupies the 
whole root interval J„. 

The existence of the separating points (^), is the key to obtaining the explicit form of 
the periodic-orbit expansion for individual roots /c„. Multiplying both sides of (|39|) by k and 
integrating from kn-\ to kn we obtain 

where we used (0) for the integral over kp and defined tUp = TrSp/^o. Since all the quantities 
on the right-hand side of (Q) are known, this formula provides an explicit representation of 
the roots k^ of the spectral equation ( ]29|) in terms of the geometric characteristics and the 
classical properties of the graph. 

In [0 a mathematical proof is presented that assures us that ( |5^ converges. In addition 



it is proved in |]I8[ that ( |5^ converges to the exact spectral eigenvalues. Both convergence 
and convergence to the correct results are illustrated with the help of a specific example in 
Sect. IV. It is also proved in |^ that the series ( |5^ is only conditionally convergent. This 



means that for proper convergence the ordering of the terms in is important. Proper 



convergence of (|54 ) is obtained if the terms in (|5^) are ordered according to the code lengths 



of the periodic orbits |18|. In other words, the sum in (p4) is to be extended over all periodic 



orbits with code lengths smaller than or equal to / which yields the approximation to 



kn- Then, on the basis of the results obtained in ||18[, we have lim^^oo ^n(0 — ^n- This 
means that (Q) is exact. It is important to note here that the ordering of terms in (|54D is 
not according to their action lengths, but according to the lengths of the code words that 
code for the periodic orbits. This is intuitively understandable, since the code length I is 
connected to the power n of the ^-matrix in (^) according to I = n/2. 

The expansion (|5^ provides an explicit representation of the roots of the spectral equa- 
tion (pPp in terms of the geometric characteristics of the graph. In a similar way one can 



obtain explicit expansions for any power of the energy levels, fc™, or any function of the 
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eigenvalues / 



IV. EXAMPLES 



The coefficients Ap in ( p4[ ) assume a particularly simple form in the case of linear graphs 
with zero sources, = 0, i = l,...,Nv- Both the vertices and the bonds of a linear 
graph can be naturally labeled by means of a single index such that Bi^2 = Bi, i?2,3 = 
B2,...,Bnv-i,Nv = B^v-i (see Fig. 2). The scaling coefficients for the momentum of the 
particle are correspondingly /5i2 = Pi, ^2,3 = P2,---,Pnv-i,Nv = Pnv-i^ ^^e bond lengths 
are L12 = Li, L2,3 = L2,...,LNy-^i,Nv = Ln^-i, the potentials are = Ui, f/2,3 = 
U2,..-,Unv-i,Nv = Unv-1 and the reduced bond actions are S^2 = = PiLi, 6*23 = 5° = 
P2L2,---,S%^_^ = S%^_-j^ = (3Ny-iL]S!y-i, respectively. In this case, if a prime periodic 
orbit p undergoes a* reflections from a vertex Vi and 27^ transmissions through it, the 
weight coefficient in the expansion (|39| ) is 

A, = X{^^-rlYK (55) 

i 

where rj is the reflection coefficient from the vertex V^, and the product is taken over all the 
vertices encountered by the orbit p. If a particle reflects from the vertex Vi traveling along 
the bond B^, the reflection coefficient is 

ri = ^'"^ , , i = 2, ...,iVv'-l, ri = -1, r^v^ = -1. (56) 
Pi-i + Pi 

We assumed Dirichlet boundary conditions at the left and right dead ends of the graph. 
The reflection coefficient changes its sign if the reflection happens from the side of the bond 
If, for a given orbit, the total number of reflections with < is Xp, then 

Ap = {-l)^^Y[\nrH'^-r^y^- (57) 

i 

The two-vertex linear graph is trivial and corresponds to a quantum particle in a square-well 
box. A quantum particle moving in a scaling step potential as shown in Fig. 3, gives rise to 
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the simplest non-trivial graph, the scaling three-vertex linear graph, shown on the bottom 
of of Fig. 3. In this case there is only one non-trivial reflection coefficient, 

r - (68) 

All the periodic orbits of the three-vertex linear graph shown in Fig. 3 correspond one to 
one with words formed from a binary code with two letters C and TZ [p"7| -p^ , where C stands 
for a reflection of the orbit off the left-most vertex (left-hand potential wall), and TZ stands 
for a reflection off the right-most vertex (right-hand wall). Thus the C, TZ code is unique 
and complete. For this system the spectral equation is 

sin(S'^A;) - sm{VLik) = 0, (59) 

where Sq = Sl + 5*2 is the total reduced action of the graph and Qi = S^— 82- With ai = r2 
and 7o = 7i = 7r/2, ( ^9]) is of the form (|4^), (^) and the number of cos-terms in (j){k) (in 
this case one term) complies with the estimate (^6]). Because of |r2| < 1, it is the spectral 
equation of a regular quantum graph. Using the explicit form (0) of the coefficients Ap in 
the expansion (0), we obtain the explicit series expansion for every root A;„ of (^9]). Thus 
the spectrum of the scaling step potential shown in Fig. 3 may be calculated explicitly and 
analytically with the help of (|S3). This by itself is a considerable advance in the theory 



of simple one- dimensional quantum systems which up to now could only be solved using 



graphical or numerical techniques 1 14, 20, 21 



We illustrate the method and the convergence of the series expansion ( ^If ) with the 
following concrete, dimensioned example of the scaling step potential of Fig. 3. Choosing 
b = 0.3, Ai = 0, and A2 = 1/2, we computed the solutions ki, kiQ and fcioo of (^) using 
three different methods: (i) exact numerical, (ii) explicit periodic-orbit expansion (|54D of the 
individual eigenvalues and (iii) numerical integration using the S-matrix representation ( |33D 
of the density of states. Addressing (i) we obtained the exact numerical values of the three 
selected roots of (|5|). The result is: ki = 4.107149, kio = 39.305209 and kioo = 394.964713. 

Turning to method (ii) we re-computed these three eigenvalues using ( |5^ directly in- 



cluding progressively more periodic orbits in the expansion (^). The result is presented in 
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Table 1 which shows the values of ki, kio and fcioo computed with (|5^) including periodic 
orbits coded by binary words of length / = 5, 10, 15 and 20, respectively. This corresponds 
to including 23, 261, 4807, and 111321 periodic orbits in the expansion (|5^, respectively. 
We observe that the accuracy does not improve monotonically, but that there is a definite 
overall improvement of accuracy with the code length. As a matter of fact, as discussed 
above and shown mathematically in [|l8l, the series ( ^4|) converges, and converges to the 
exact results of /c„ in the limit of / — oo. 

Turning to method (iii) we note that due to the exponential proliferation of periodic 
orbits, it becomes progressively more difficult to compute the codes of longer periodic orbits. 
Nevertheless, with the help of a numerical procedure, we are able to illustrate the convergence 
behavior of (|5^) for code lengths much longer than / = 20. Starting from (^) we compute the 
S'-matrix numerically and perform all the steps leading up to (|54D numerically. In particular 
this method involves numerical computation of S'-matrix powers and numerical integration 
over k. Within any given level of numerical accuracy this method is completely equivalent 
to the method of summing the orbits, but allows us to extend the computations such that 
we effectively include all periodic orbits up to code length / = 150. This corresponds roughly 
to 2^^*^ = 1.4 X 10'^^ periodic orbits, since the periodic orbits on the three- vertex linear graph 
are coded by a binary code. This estimate is substantiated by an analytical estimate of 
the number of periodic orbits. For the three-vertex graph the periodic orbits are binary 



necklaces over the two symbols C and TZ [jT8| . The number of binary necklaces of length i 
is given by p2 

m) = 7 E <f>in) 2'/^ (60) 

where the symbol "n|£" denotes "n is a divisor of f , and 0(n) is Euler's totient function 
defined as the number of positive integers smaller than n and relatively prime to n with 



0(1) = 1 as a useful convention. An upper limit for A is obtained if we use (pO| ) in the case 
where £ is a prime number. In this case (pO) reduces to 



Af{p) = i [0(1) 2^ + 0(p) 2^] = ^ [2^ + 2(p - 1)] , (61) 
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where p is prime. Thus, in the hmit of p — oo we have A^(p) — > and therefore, as an 
upper hmit 

A=.tol« < hml!«l = ,„(2). (62) 

p— >oo p p^oo p 

Thus, according to this estimate, the total number of periodic orbits of length 150 is again 

#(150) ~ e^5°^ = 2i5°. (63) 

We also computed numerical estimates of A. Using the exact formula (^) for counting 
periodic orbits in (|^) and including periodic orbits with code lengths of up to / = 1000, we 
found A > 1.987, consistent with the estimate (l6^) . For I = 150, relevant for our numerical 
example, the asymptotic regime is not yet reached and we find A ^ ln(1.943). This value for 
A can be used for a more refined estimate of the number of periodic orbits of length / = 150, 
#(150) ~ 1.943^^° ^ 2 X 10^'^. Clearly, computing the codes of that many periodic orbits 
and summing them up in (|5^) is beyond the storage capacity and power of any currently 
existing computer, but is apparently no obstacle to the numerical simulation of that many 
periodic orbits included in (0). Figure 4 illustrates the rate of convergence of the eigenvalues 
ki, kiQ and fcioo obtained with method (iii) as a function of code length / up to / = 150. 
Shown is the relative error e; = |A;„(/) — kn\/kn for = 1, 10 and 100. The error is seen 
to decrease on average as a function of increasing periodic-orbit length /. From Fig. 4 we 
obtain approximately e; = \kn{l) — kn\/kn ~ 1//^ on average. 
For a four-vertex linear graph, the spectral equation is 

sin(S'Q/c) = r3 sin(f2iA;) — r2r3 sm{fl2k) + ^2 sin(n3A;), (64) 

where = 5° + + Si Qi = + - Si Q2 = S^ - S^ + ^3°, = 5? - S^ - S^ and 
r2, r3 are the reflection coefficients at the vertices V2 and V3, respectively. With 7j = 7i/2, 



i = 0,1, 2, 3, this spectral equation is of the form (^4|), (^) and the number of cos-terms in 
(f){k) (three in this case) complies with (|46|). For 



ksl + 1^2^31 + ^2! < 1, (65) 
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the four-vertex linear graph is regular. In this case the energy values of the four-vertex 
linear graph may be calculated exactly using the periodic-orbit expansion (p^. According 
to Fig. 5 the set of r2, values that fulfill ( |65D occupies a diamond-shaped area bounded 
by the functions rs = ±(1 — |r2|)/(l + \r2\)- This observation proves that regular quantum 
graphs are an important, finite- measure subset of quantum graphs. 

The set of regular quantum graphs is much wider than the three- and four- vertex quantum 
graphs discussed in detail above. Since, as indicated by ( ^9|) and (|64|), the amplitudes 
Qi in (^51) involve products of vertex reflection coefficients, and since the vertex reflection 
coefficients of a linear quantum graph (via proper choice of the bond potentials) are free 
parameters of the quantum graph, a finite-measure set of regular quantum graphs exists for 
any given linear graph. It is possible that more complex graph topologies, such as rings 
and stars, may also admit a set of regular quantum graphs. This topic is currently under 
investigation. 



V. DISCUSSION, SUMMARY AND CONCLUSION 

In this paper we defined and studied a subset of quantum graphs: regular quantum 
graphs. Regular quantum graphs satisfy the regularity condition (^) which implies that 
the roots of the spectral equation are confined to regularly spaced root intervals. One 
and only one root is found per root interval. This property allows us to derive rigorous, 
converging periodic-orbit expansions for individual energy levels of regular quantum graphs. 

We often hear the comment that the expansion ( ^If ) for kn cannot possibly converge, 
since (p9D is divergent. This comment is invalid. There is a fundamental difference between 
(^) and (|54D. Equation (p9D is a periodic-orbit expansion of the kernel of a functional (a 



series of Dirac delta "functions"), whereas ( |54D is a periodic-orbit expansion of a simple c- 
number. On the level of (^) the concepts of convergence or divergence are not even defined. 
Only after multiplying (^) with a test function (a function with compact support |23|) and 
integrating over k, are the concepts of divergence and convergence defined. In this sense 
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even (^) is convergent. This is also known as convergence in the distribution sense p3[ and 



leads to proper convergence in the usual sense of elementary undergraduate-level analysis 



after multiplying with the test function and integrating. Thus ([5^ ) converges in the usual 
sense of elementary calculus since our method of multiplying ( |5^ ) with k and integrating 
from kn-i to kn is equivalent with multiplying (BSI) with the test function 



1, for A;„_i < k < kn, 
g{k) = <{ ' « 1 - - (66) 

0, otherwise, 



and integrating over k. Since g{k) is a proper test function the convergence of (|54D is 
no longer a mystery. 

Another more serious comment concerns the sense in which quantum graphs are clas- 
sically chaotic. Quantum graphs are based on a one-dimensional network of vertices and 
bonds. Therefore a dynamical Liapunov exponent [0,0] cannot be defined in the traditional 



sense of exponentially diverging initially close trajectories. As explained in 0, however, this 
is no obstacle to associate a classical phase space with a quantum graph and to show that 
the classical dynamics in this phase space is mixing [Q]. Moreover, quantum graphs fulfill 
another property of quantum chaos, the exponential proliferation of classical periodic orbits 
as manifested by a positive topological entropy (see Sects. Ill, IV). Therefore, quantum 
graphs have been called "paradigms of quantum chaos" 0|. Inasmuch as the positive topo- 
logical entropy is concerned regular quantum graphs qualify as quantum chaotic systems. 
Regular quantum graphs do not show all the characteristics of "fully developed" quantum 
chaos. For instance, because of the existence of the forbidden zones i?„ (see Sect. Ill) they 
definitely do not show a Wignerian nearest neighbor statistics. We do not believe that this 
is a problem since there is no universally accepted rigorous definition of quantum chaos 
that requires "Wignerian statistics" as one of the necessary conditions. The only broadly 
accepted criterion is that "quantum chaos" deals with quantum systems that are chaotic 
in their classical limit. Based on this criterion, together with the mixing property and 
the positive topological entropy (see Sect. IV) regular quantum graphs definitely qualify as 
quantum chaotic systems, i.e. systems chaotic in the classical limit. 
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A more delicate point concerns the derivation of ( ^41) by integrating the fluctuating 
part in (^) term-by-term. Since the resulting periodic-orbit expansion is only conditionally 
convergent, this is reason for concern. We addressed this point from a rigorous mathematical 



point of view in |T8[. We were able to prove rigorously that the interchange of integration 
and summation in the k integral of kp{k) is allowed. Thus, term-by-term integration of 
is justified, validating the final result (|5^) . 

To our knowledge this is the first time that the energy levels of a class of classically 
chaotic systems are expressed one by one with the help of convergent periodic-orbit expan- 
sions. Studying specific examples of regular quantum graphs we proved that the class of 
regular quantum graphs is not empty; that it is in fact an important finite-measure subset 
of quantum graphs. The explicit formulae of individual quantum energy levels obtained in 
this paper remind us of the Einstein-Brillouin-Keller (EBK) method for the quantiza- 
tion of integrable classical systems. But there are important differences. Regular quantum 
graphs do not correspond to classically integrable systems. In fact, due to the importance 
of non-Newtonian periodic orbits, the number of classical periodic orbits proliferates expo- 
nentially with the code-length. This proves that even regular quantum graphs, as defined 
by the regularity condition (^7\), are classically chaotic systems with a positive topological 
entropy. For one of our examples, the scaling step potential discussed in Sect. IV, we com- 
puted estimates of the topological entropy analytically and numerically. In both cases it 
turned out to be positive and close to ln(2). This proves that, at least for the cases studied, 
the classical limit is chaotic. Another difference to EBK theory is that the periodic orbits 
in regular quantum graphs are not confined to phase-space tori. Finally, in contrast with 
EBK theory, a semiclassical theory which does not usually return exact results, our formulae 
are mathematically exact. In summary, despite the apparent complexity and exponential 
proliferation of the periodic orbits of regular quantum graphs, the organization of the roots 
of the spectral equation into regularly spaced intervals makes it possible to pinpoint every 
single energy eigenvalue of a regular quantum graph analytically and exactly by an explicit, 
convergent periodic-orbit expansion. 
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Figure Captions 
Fig. 1: A generic non-planar graph with six vertices and ten bonds. 

Fig. 2: An example of a (Manhattan) step potential (a) and its associated linear graph (b). 
Also shown is a non-Newtonian periodic orbit characterized by six above-barrier reflections. 

Fig. 3: A scaling step potential (top), equivalent to a three- vertex hnear graph (bottom), 
as an example of a regular quantum graph. A Newtonian {CIV) and two non-Newtonian (£, 
CTZTZ) periodic orbits are also shown together with their C — TZ codes. 

Fig. 4: The deviation e; = \kn{l) — kn\/kn of the exact eigenvalues for ki, kio, and kioo from 
the corresponding values obtained via the series representation, as a function of the lengths 
/ of the periodic orbits. 

Fig. 5: Parameter space (r2,r3) of a four-vertex linear quantum graph. Parameter com- 
binations in the shaded region correspond to regular quantum graphs. This demonstrates 
that the subset of regular quantum graphs within the set of all four- vertex linear quantum 
graphs is non-empty and in fact of finite measure. 
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TABLES 



root 


1 = 5 


I = 10 


I = 15 


I = 20 


exact 


error 


ki 


4.11608 


4.11653 


4.10721 


4.10513 


4.10715 


0.00202 


kio 


39.28658 


39.29807 


39.30730 


39.30521 


39.30521 


0.00000 


hoo 


394.94770 


394.95647 


394.96622 


394.96456 


394.96471 


0.00016 



Table Captions 

Table 1: Successive approximations of the eigenvalues ki (first row), k^o (second row) and 
kioo (third row) of a specific scaling step-potential (see text for details) as a function of code 
length I (columns 2-5). The exact values of ki, kio and fcioo are listed in column 6. Column 
7 lists the absolute errors \kn{l = 20) — A;„| for n = 1, 10 and 100. 
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